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ON ' Abstract 



We exploit the fact that within massless perturbative QCD the same Green's 
function determines the hadronic contribution to the r decay width and the mo 



' ments of the e + e~ cross section. This allows to obtain relations between physical 

observables in the two processes up to an unprecedented high order of perturbative 
QCD. An precision measurement of the r decay width allows one then to predict 
the first few moments of the spectral density in e + e~ annihilations integrated up 
to s ~ m 2 . with high accuracy. The proposed tests are in the reach of present 
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experimental capabilities. 



Due to the availability of ever-increasing orders of perturbative QCD results on the 
one hand and continuing improvements on the quality of experimental data on the other 
hand, high order precision test of the QCD part of the Standard Model are now becoming 
feasible [11] . In contrast to the electroweak sector where the expansion parameter for the 
perturbative calculations is sufficiently small and one or two corrections provide a theoret- 
ical accuracy sufficient for a comparison with present experiments, the situation with the 
strong interactions is different, maybe even in principle. The relevant expansion parame- 
ter is large, such that a desired accuracy of order 1% requires four- to five-loop accuracy 
in the perturbative expansion. This is close to the border of today's computational pos- 
sibilities. In addition, with such a large number of terms, one may already encounter the 
asymptotic nature of the perturbation series in which case no further increase of precision 
is possible in straightforward perturbation theory. 

The above question of numerical convergence is influenced to a large extent by the 
renormalization scheme dependence of the truncated perturbation series |2|, || ||. There- 
fore it is desirable to obtain predictions for observables which are renormalization scheme 
independent ||. 

In this letter we propose a new test of perturbative QCD that is free from renormaliza- 
tion scheme ambiguities and is of higher precision in terms of the number of corrections 
than any previous test. This test directly expresses one observable through another ob- 
servable using massless perturbative QCD as the underlying theory. A disagreement with 
experimental data would imply a failure of the applicability of perturbative QCD. 

The advocated comparison is made for the moments of the spectral density in e + e~ 
annihilation and the hadronic contributions to T(r — > v T + hadrons). Within massless 
perturbative QCD these quantities have been computed up to four-loop order. We define 
a reduced semileptonic r decay ratio r T through 

r( p^ + hadxons) 
r(r -»• v T + 11 + Vp) 

where R® is the lowest order partonic value of R T . To lowest order in the electroweak 
interactions and assuming \V U( i\ 2 + \ V US \ 2 = 1 it takes the value R® = N c = 3. Corrections 
due to higher order electroweak interactions and power suppressed terms are small || [7]] 
and can be easily accounted for by rescaling the partonic ratio R^. Thus the reduced 
decay width r T is determined by massless perturbative QCD, for which the axial and 
vector contributions are identical. The expansion for r T starts directly with a(/i 2 ), where 
a = a s /ix. The dependence of a(/i 2 ) on the renormalization scale /i is determined by the 
renormalization group equation 

/U 2 T^T = P( a ) = -° 2 (A) + Pia + P 2 a 2 + p 3 a 3 + ...). (2) 
In minimal subtraction schemes one has 

A) = ^ = 2.25, 

* = iir 4 ' 

(3 2 = ^^ 10.060, 
H 384 

140599 445 w , 

03 = + C 3 w 47.228 3 

M3 4608 32 sv ; v ; 
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for N c = rif = 3. Here the recently computed four-loop coefficient (3$ || is also included. 
In e + e~ annihilation the cross section is determined by the imaginary part of the vacuum 
polarization, 

R e+e -(s) = 127rImn( S ) = N C J2QK^ + r(a)) = 2(1 + r(s)). (4) 
In perturbative QCD one has 

r(s) = a(fi 2 ) + (h + (3 L)a 2 (fi 2 ) + (k 2 - fvr 2 /? 2 + (2fah + fa)L + /3 2 L 2 )a 3 (/i 2 ) 
+ (A; 3 - tt 2 /^ - §7r 2 /3 /3i + (3/3 A; 2 + 2fak x + fa - tt 2 $)L 
+(3(3 h + |A)L 2 + /5 3 L 3 )a 4 (/i 2 ) + . . . (5) 
with L = ln(/i 2 / s) and 

h = —- 9C(3) * 1.63982, 

*. - ^Z-^C(3) + fc(5)« 6.37101 ,6) 

in the modified minimal subtraction (MS) scheme |J, while k$ is still unknown. 
We can define moments of the reduced part r(s) of the spectral density, 

r B (« ) = (n + l) [ S °-(-) n r(s) (7) 



/o s \s , 

which, for the sake of convenience, are normalized to unity for r(s) = 1. In terms of these 
moments the reduced decay width ratio r T is given by 

r T = 2r (m 2 ) - 2r 2 (m 2 ) + r 3 (m 2 ). (8) 

Eq. (|8[) can be inverted within perturbation theory. One can then express the perturbative 
representation of one observable, i.e. any given e + e~ moment r n (m 2 ), in powers of r T using 
the perturbative expansion of the r decay observable. The strong coupling constant a s 
in any given scheme (MS for the present letter) serves only as an intermediate agent 
to obtain relations between physical observables. The reexpression of one perturbative 
observable through another is a perfectly legitimate procedure in perturbation theory and 
is standard practice in electroweak theory. The result is independent of the choice of the 
renormalization scheme. Solving for the moments r n (m 2 ) in terms of the reduced r decay 
width r T one finds 

r„(m 2 ) = f 0n r T + f ln r 2 T + f 2n r 3 T + f 3n r 4 T + f 4n r 5 T + 0(r 6 T ), (9) 



where 




fon 


= /(0,n), 


fin 


= A)I(l,n), 


fin 


= /3 2 (/(2,n)+p 1 J(l,n)), 


f3n 


= /? 3 (/(3,n) + (J T (2) - J T (1) 2 - \<K 2 )I{l,n) + 


fin 


= /? 4 (/(4,n) - 3(J T (2) - / r (l) 2 - |7r 2 )/(2,n) 




+2(/ T (3) - 3/ T (l)J T (2) + 2J r (l) 3 )J(l, n 




+pi(f 7(3, n) + 5(/ T (2) - J r (l) 2 - i7r 2 )/(l, : 
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\pj(2, n) + p 2 1(1, n 



(10) 



with 

ml 

(n + 1) TO ' 

2I(m, 0) -2I(m,2)+I(m,3), 
m (m\ 

I(m) + J2[ }l T (p)I(m-p,n). (11) 
p =o \ P / 

The coefficients J(m, n) used in Eqs. (|T0|) can be determined iteratively by solving Eq. (|il"|) . 
In general, the coefficient functions /j n depend on so/m 2 , but here we limit our analysis 
to the value s = m 2 T . The pj are scheme independent quantities given by 

pi = ^ 0.79012, 

Po 

92 = 4 - A^i + /3o(k 2 ~ kl)\ « 1.03463, 
Po 

p 3 = [/? 3 -2/5 2 A; 1 +/? 1 fc 1 2 + 2/?o(A;3- 3^^2 + 2^)] « 0.17558fc 3 - 2.9795. (12) 

Po 

The coefficient /4 n contains the unknown five-loop coefficient fc 3 . Later on we present 
an estimate of the coefficient k 3 in order to be able to assess the importance of the r\ 
contributions to the perturbation series. 

We now turn to the numerical analysis. From the experimental value R° xp = 3.649 ± 
0.014 (as a recent review see e.g. J7|) for the semileptonic r decay ratio one gets r^ xp = 
0.216 ± 0.005. This value is obtained from Eq. ([!]) using R° T = 3. Nonperturbative and 
electroweak corrections lie within the error bars. In order to investigate the convergence 
properties of the series we calculate the first few moments by inserting the numerical 
values of all known coefficients in Eq. @. One has 



r-1/2 


= 0.216(1 + 0.203 


+ 0.786 + 2.637 + (10.601 + 0.0040fc 3 )) + 0(r°), 


(13) 


r 


= 0.216(1 


- 0.284 


- 0.069 + 0.110 + (0.211 


- 0.0057fc 3 )) + 0(4), 


(14) 


ri 


= 0.216(1 


-0.527 


- 0.143 + 0.177+ (0.361 


-0.0106*3)) +0(r% 


(15) 


r 2 


= 0.216(1 


- 0.608 


- 0.115 + 0.269 + (0.433 


-0M22k 3 )) + O(r 6 T ), 


(16) 


r 3 


= 0.216(1 


- 0.648 


- 0.091 + 0.317 +(0.444 


-0M31k 3 )) + O(r 6 T ). 


(17) 



We have started with the moment r_i/ 2 because it represents the cross section integrated 
over the energy y^s which would be the natural choice of moment when considering the 
usual representation of the experimental data. It is quite apparent from Eq. ([i3|) that 
the series for n = —1/2 does not converge at all. A possible explanation is that the low 
energy region of the spectral function is emphasized too much for perturbation theory to 
work properly. 

To discuss the convergence properties of the perturbation series for the other moments 
(n = 0,1,2,3) we distinguish two scenarios depending on how the unknown r\ term is 
accounted for. These are 

(a) Nothing is known about the r\ term (i.e. about k 3 ). Then we truncate the series at 
the smallest term and consider the next term as an estimate of the error. We call 
this the conservative estimate. 



I(m,n) = 
I T (m) = 
I(m,n) = 
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(a) 


(b) 


(c) 


Rq 


2.28 ±0.05 


2.30 ±0.02 


2.15 


Ri 


2.14 ±0.08 


2.18 ±0.04 


2.06 


Ri 


2.12 ±0.12 


2.18 ±0.06 


2.00 


Rs 


2.11 ±0.14 


2.18 ±0.07 


1.99 



Table 1: Numerical values of the e + e moments for the conservative estimate (a), for the 
educated estimate (b), and for the experimental e + e~ moments (c). 



(b) We adopt a value k 3 = k\jk\ ~ 25 as estimated by using Pade techniques [|TI], 
|i~lf . There are other methods to estimate k 3 which give results close to the Pade 
estimate |T2[ . With this choice of k 3 the series appears to consist of two alternating 



series, one of them with large terms and the other with small terms. With such a 
structure our best prediction lies between the third order and the fourth order partial 
sum, so the error is half of the difference of the two. Note that the correction to the 
fourth order partial sum is always small. We call this the educated estimate. 

We collect our numerical results for the full moments R n = 2(1 ± r n ) and their errors in 
Table [I| Here we also list results for the experimental moments using the experimental 
e + e~ data from a recent compilation fi~3|| . On average the error bars of the conservative 
estimate ( factor two bigger than those of the educated estimate (b). For R the 

error is very small (even for the conservative estimate) because the perturbation series 
is reliable and converges well numerically. The data seem to be systematically lower by 
about 7%. The resulting statistical error for the moments of the data is negligible, the 
systematic error could be as large as 10% [13|, For the higher moments the theoretical 
predictions are less accurate. The experimental values lie within the error bars of the 
conservative estimate (a). If one takes the educated estimate (b), the experimental values 
are consistently lower than the theoretical prediction. If the data points were increased 
by 7%, the moments would be consistent with both theoretical predictions. 

We have also checked on the non-perturbative contributions to the moments Rq, Ri 
and i?2- The non-perturbative contribution to Ro vanishes, since there is no gauge invari- 
ant operator of dimension two. The non-perturbative contributions to R\ and R2 are well 
within the truncation error of the perturbation series. 

As a last point we consider the relation in Eq. which holds to arbitrary order in 
perturbative QCD. If we use experimental data on both sides of Eq. (||), we find the value 
3.449 for the right hand side, so the compiled data of the e + e~ annihilation are lower by 
6% in comparison with r T . On the other hand, if we use the predicted moments of Table |1] 
to evaluate the right hand side of Eq. (|8|), we obtain 3.648 in both cases (a) and (b). This 
is very close to the input value r T = 3.649 ± 0.005 and within our truncation errors which 
means that our procedure is self-consistent. This observation confirms the basic result of 
perturbative QCD that moments of vector and axial vector spectral functions are equal to 
one another. If against all expectations future precise e + e~ data does not satisfy Eq. (||) 
then this would imply that the moments of the vector and axial vector spectral functions 
differ and that perturbative QCD is in trouble at the scale of w? r . 

In summary we have established perturbative relations between different sets of ob- 
servables with one more term in the perturbation series than ever considered before. This 
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seemingly surprising result can be traced to the fact that the observables are generated 
from the same Green's function. The coefficients fi n in our perturbative expansion are 
scheme independent quantities. The stated relations between the observables are quite 
useful and can be tested since the observables can be measured directly, independently 
and very precisely. 

Our analysis could be extended in two directions as concerns possible improvements 
on the convergence properties of the series expansion. One can either vary the endpoint sq 
of the moment integration in order to optimize the series expansion in Eq. (£5|) or one can 
attempt to find an optimal linear combination of moments that has improved convergence 
properties without a concomitant loss of precision through cancellation effects. 

In view of the fact that the e + e~ cross section will be remeasured with considerably 
improved accuracy at DA$NE we suggest the following treatment of the data. For 
our analysis only integrals of the data are needed. Excluding the resonance region (say 
E < 1 GeV) the resulting spectral density is sufficiently smooth because the integral is a 
linear combination of data points. It is therefore sufficient to have only a few points but 
measured with very high precision. 
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